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Abstract 

We calculate the processes e + e~ — > Atc and r — > z/ T 47r to 0(p 4 ) in the low-energy expansion 
of the standard model. The chiral amplitudes of 0(p A ) can be extended via resonance 
exchange to energies around 1 GeV. Higher-order effects have been included in the form 
of u, ai and double p exchange and by performing a resummation of the pion form 
factor. The predicted cross sections and the branching ratios BR(p° — > An) are in good 
agreement with available data. 
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1 Introduction 



Electron positron annihilation into hadrons has played an important role in the develop- 
ment of modern particle physics. Precise knowledge of the cross section is essential for 
many purposes, in particular for the determination of the hadronic contribution to the 
anomalous magnetic moment of the muon and for running the fine-structure constant up 
to Mz to analyse electroweak precision measurements. 

At high energies, the inclusive cross section can be calculated in QCD. It provides 
one of the standard tests of QCD allowing for the extraction of the strong coupling 
constant. At energies below approximately 2 GeV, the different exclusive channels are 
measured separately. As perturbative QCD cannot be applied to those exclusive processes 
the theoretical challenge consists in modeling them in a way that is at least consistent 
with QCD. 

At very low energies (E << 1 GeV), the most reliable approach is furnished by chiral 
perturbation theory (CHPT), the systematic low-energy effective theory [1], ^ |3| of the 
standard model. Although the low-energy expansion of CHPT breaks down at typical 
hadronic scales of 0(M P ) it can still provide important constraints how to match the 
low-energy amplitudes on to the intermediate-energy region governed by meson resonance 
exchange 0, 0], A simple but very illustrative example is the pion form factor measured 
in e + e~ — > 7r + 7r~ that can be continued from threshold to the region beyond 1 GeV in a 
straightforward way |], [7], §] . 

Inspired by this success and by the phenomenological importance of four-pion pro- 
duction, especially for the calculation of a(E), we have undertaken a systematic study of 
e + e~ — > 4ir in CHPT (with the two possible charge configurations 27r°7r + 7r~, 27r + 27r~). At 
first sight, this does not appear very promising with the threshold E = AM^ ~ 560 MeV 
already in the vicinity of the most prominent meson resonance, the p meson. However, 
once again the chiral amplitudes of 0(p A ) can be continued into the resonance region 
and the decay rates T(p — > 4n) can be calculated with reasonable accuracy. We will also 
consider the energy dependence of cross sections for E < 1 GeV. 

There is an essential difference between the two- and four-pion modes. Whereas the 
two-pion amplitude is completely dominated by p exchange even beyond 1 GeV the 
situation is much more involved for four-pion final states where at least p, uj and a\ 
exchange are relevant (e.g., Refs. |9], [ID, |TT|]). It is then all the more important to have 



unambiguous theoretical guidelines for the construction of those amplitudes such as the 
correct low-energy behaviour dictated by QCD. 

In addition to electron positron annihilation, multi-pion final states can also be studied 
in t decays (for reviews of the theory see, e.g., Refs. [|T2|, [13]] )• In the limit of isospin sym- 
metry that will be assumed throughout this paper both the two- and the four-pion modes 
are related. There is again an important difference between the two modes. Whereas the 
annihilation amplitude and the decay amplitude for two pions in the final states are in 
one-to-one correspondence the situation is more subtle in the four-pion case [13[|: given 
the amplitude for e + e _ — > 2tt°7t + 7t~ or for t~ — > u T 27i^7T + 7i , the three remaining an- 
nihilation and decay amplitudes are uniquely determined but not vice versa. Therefore, 
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in the isospin limit it is sufficient for a complete determination of all four amplitudes to 
construct the amplitude for the 27r°7r + 7r~ channel only. All calculations in this paper will 
be performed for this particular channel. 

In Sec. |2| we collect the kinematics, matrix elements and cross sections for the process 
e + e~ — > Atc. We recall the isospin relations relating e + e~ annihilation and r decays into 
four pions. To make Bose symmetry and C invariance manifest, we express all matrix 
elements in terms of a reduced amplitude that reduces the size of the various amplitudes 
roughly by a factor four. The leading-order amplitudes of 0(p 2 ) for both e + e~ annihila- 
tion and t decays are presented in Sec. [3[ The matrix elements of 0(p 4 ), consisting of 
both one-loop and tree-level contributions, are calculated in Sec. |j. The structure of the 
local amplitude of 0(p 4 ) determines the resonance exchange amplitudes generated by p 
and scalar exchange. The relevant chiral resonance Lagrangian and the resulting matrix 
elements are presented in Sec. [5| To extend the amplitudes into the resonance region, 
additional contributions are necessary. In Sec. ^| we analyse double p, to and a\ exchange 
to obtain the final amplitude. We collect the leading terms in the low-energy expansion 
of the pion form factor and we resum those terms to the complete p dominated form 
factor. The energy dependent cross sections for the two channels and the partial widths 
r(p° — > 27r°7r + 7r~), r(p° — > 27T + 27r _ ) are analysed in Sec. [7[ We compare our results with 
available data in the region 0.65 < E'(GeV) < 1.05 and we discuss the necessary steps for 
proceeding to higher energies. Our conclusions are summarized in Sec. Three appen- 
dices contain a discussion of the isospin relations, a brief summary of the two possibilities 
for incorporating spin-1 mesons in chiral Lagrangians and a collection of numerical input 
for the calculation of cross sections. 

2 Kinematics and symmetries 

The amplitude for the process 

e + (k + )e~(k_) -> 7r(pi)7r(p 2 )7r(P3)7r(P4) 

is written in the form 

e 2 

M = . v(k + )<y t ,u(k_)J fM (p u p 2 ,p 3 ,P4) , q = k + + k_ 
q z + ie 

with J M the pionic matrix element of the electromagnetic current 

J^(Pi,P2,P3,Pa) = (7r(Pi)7r(P2)vr(p 3 )7r(p4)|^dm( )l ) 
The differential cross section is then given by (setting m e = 0) 

do = -^rsifl ^)5 (4) (g -YpiY v JJI 



Eft , (2-1) 



(2.2) 



(2.3) 
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with the leptonic tensor 

a 2 

= h%k v _ + ktk\ - —g^ . (2.4) 

We will only be interested in the integrated cross sections cr(g 2 ) where the appropriate 
statistical factors have to be applied for the two possible channels 27r°7r + 7r~ and 27r + 27r~. 
With the charge assignments 

rc (p 1 )TC (p 2 )-iT + (p 3 )TT-(p 4 ) , 7r + (p 1 )7r + (p 2 )7r"(p 3 )7r"(p4) , 

a convenient set of Dalitz variables^ is 

Q 2 , s = (px+ p 2 f , v = (p 3 - p A ) ■ (p! - p 2 )/2 , 

U=Pi-q (i = 1,...,4) . (2.5) 

4 

There is a redundancy due to the relation £j = q 2 but especially for displaying symme- 

tries of the amplitudes it is useful to keep the complete set. For compactness of notation, 
we will often express amplitudes in terms of the various scalar products instead of using 
s and v. 

In the isospin limit that will be assumed throughout this paper, the current matrix 
element for the 27r + 27r~ channel can be expressed in terms of the matrix element for the 
27r°7r + 7r _ channel Hl3|: 



(vr (pi)7r (p 2 )7r + (p 3 )7r-(p4)|J c 1 m (0)|0) := J"( P i,p 2 ,P3,p 4 ) ( 2 -6) 
(7r + (pi)7r + (p 2 )vr _ (p3)7r _ (p4)|J c 1 m (0)|0) = J M (pi,p 3 ,P2,p 4 ) + J^(pi,P4,P2,Ps) 

iP3,P1,Pa) + J^{P2iP4:iPl-,P-i) ■ (2.7) 

Likewise, the matrix elements of the charged vector current relevant for r decay can also 
be expressed in terms of J fJ, (pi,P2,P3,P4] 



(n-(p 1 )n-(p 2 )n + (p + )n°(p )\V c ^0)\0) = V2 {J^p + , Pl ,p 2 ,p ) + J^p + ,p 2 , Pl ,p )}(2.8) 
(7r (p 1 )7r (p 2 )7r (p 3 )7r-(p-)|y c ^(0)|0) = V2 {J^ Pl ,p 2 ,p.,p 3 ) 

+ J"iPl,P3,P-,P2) + J'"{p2,P3,P-,Pl)} , (2.9) 

with the usual normalization of the charged vector current V£ = dj^u in terms of quark 
fields. We come back to these matrix elements in Sec. [| for the lowest-order chiral expan- 
sion. 

In addition to isospin, the electromagnetic current matrix elements are also con- 
strained by gauge invariance, Bose symmetry and C invariance. It is sufficient to con- 
sider these symmetries in the 27r°7r + 7r _ channel. Via the isospin relation (|2.7|) , the matrix 
element for the 27r + 27r~ final state is then automatically gauge invariant, Bose symmet- 
ric and odd under C. Of course, all these symmetries are implemented in CHPT so the 



lr rhe main convenience is in making symmetries manifest as discussed below. 
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Figure 1: Tree diagrams for 7* — > Arc. In this and subsequent figures solid lines denote 
pions; the wavy line stands for a virtual photon. 

following relations emerge automatically in the calculation and need not be imposed a 
posteriori. 

Gauge invariance (vector current conservation) implies 

q^(Pl,P2,P3,P4)=0 . (2.10) 

The remaining symmetry constraints for J^{pi-,P2-,P3iPi)i 

Bose symmetry: J fl (Pi,P2,P3,Pi) = J^{P2,Pi,Pz,P4) (2.11) 
C invariance: J t *(pi,P2,P3,P4) = -J (l (pi,P2,Pi,P3) , (2.12) 

can always be made manifest by writing 

J"(P1,P2,PS,P4) = A fl (p 1 ,p 2 ,P3,P4) + A' 1 (p2,Pl,P3,P4) 

-A»{p 1 , P2 ,Pi,p 3 ) - A' i (p 2 ,pi,p4,p 3 ) (2.13) 

in terms of a reduced amplitude A fJ, (p 1 ,p 2 ,P3,P4) that is not further constrained by C 
invariance or Bose symmetry. In this paper, we shall express all amplitudes in terms of 
A fl (pi,p2,p3,p4 1 ). This makes the sometimes quite elaborate matrix elements considerably 
more compact. Another simplification consists in dropping terms in ^4 M (pi,P2,P3,P4) and 
therefore also in J fJl {pi,P2,P3,P4) that are proportional to g M . Of course, such terms cannot 
contribute to the differential cross section ( |2.3|) . Dropping such terms may lead to seeming 
violations of gauge invariance. It goes without saying that those terms can always be 
recovered uniquely for a given matrix element by imposing current conservation (|2.10|) . 
This trivial remark will be relevant when calculating r decay matrix elements via the 
isospin relations fl2.8y2.9|) . 



3 Amplitudes at leading order 

At leading order in the low-energy expansion of the standard model, 0(p 2 ), the amplitudes 
for e + e~ — ► 4n are determined by "virtual" bremsstrahlung. The corresponding diagrams 
shown in Fig. [I] are easily calculated from the chiral Lagrangian of 0(p 2 ) for chiral SU (2) 

i 

£2 = ^(DpUDtrf + X U ] + X ] U) . (3.1) 
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The notation is standard (see, e.g., Ref. [0). For our purposes, the covariant derivative 
of the pion matrix field U contains only the external electromagnetic field . The scalar 
field x is proportional to the light quark mass rh (we set m u = := rh) and (...) 
denotes the 2-dimensional trace: 

D»U = <W + ~eA*[T 3 ,U} 

X = 2Bm 1 = M 2 1 = M 2 [l + 0{m)} 1 
F n = F[l + 0(m)] = 92.4 MeV 
(0\uu\0) = -F 2 B[1 +0(m)} . (3.2) 

As discussed in Sec. ||, we express our results in terms of the basic amplitude 
A^(pi,P2,P3,Pa) defined in ( |2.13| ) that determines all matrix elements of interest (|2.6|). 
. . . , ( p.9|) . For the tree-level amplitude of 0{p 2 ) corresponding to the diagrams of Fig. [I] 
one finds 

s ~ M 2 pj 
F 2 2h - q 2 

The complete current matrix element (|2.6|) at 0(p 2 ) is therefore given by 



Af 2) (pi,p 2 ,P3,p 4 ) = ^ 77 ^ F3 ~, ■ (3.3) 



J &(^,P3,p 4 ) = -pj- -^rV 2 ) ■ (3 - 4) 

We have used the physical quantities M n , F v in these matrix elements. The renormaliza- 
tion of M, F to M v , F^ is an effect of at least 0(p 4 ) and will of course be included in the 
amplitudes of next-to-leading order. 

The matrix element ( |3.4|) has an obvious interpretation: (s — M 2 )/F 2 is the leading- 
order amplitude for 7t°tt — > 7t + tt~ and the second factor reduces to the usual 
bremsstrahlung factor for real photons (q 2 — > 0). Although we do not discuss r — > 4n 
decays in any detail here, we also display the tree-level current matrix elements fl2.8|) , 
(p.9|), after repairing gauge invariance in ( |3.4| ) by adding the appropriate amplitude pro- 
portional to q^: 

(vr-(pi)7r-(p 2 )7r + (p + )7r°(p )|^(0)|0) = 

V2 i -2(p» + +p%)+2R^p )p + -(q-p ) + j2R^p l ){2p -(q-p l )-M 2 }) (3.5) 



(7r ( Pl )7r Q (p 2 )7r (p3)7r-(p-)|^(0)|0) 
y/2 



8=1 



^ |4j£ - M, 2 fl"(p_) - Y.R»(p,)l2p_ ■ (q- Pl ) - M 2 W }\ (3.6) 



q 2 — 2p ■ q 



In the chiral limit (M n = 0), these matrix elements have the same structure as in 
the standard reference on the subject |TjJ. There are two misprints in Ref. |15| that have 
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Figure 2: One-loop diagrams for 7* — > 47r. The virtual photon is to be appended on all 
possible lines and vertices. Wave function renormalization diagrams are not shown. 

propagated into some of the subsequent literature: Eqs. (2.9a) and (2.9b) of Ref. |L5| must 
be multiplied by the same factor v / 2/(3v / 3)- We have checked the matrix elements (|3.5|) , 
( p.6| ) also by direct computation from the chiral Lagrangian ( |3.1| ) (adding the appropriate 
external charged vector field). 

The amplitudes of 0(p 2 ) define the low-energy limit that all amplitudes must satisfy 
in order to be consistent with QCD. By themselves, they cannot be expected to provide 
a realistic approximation to the physical amplitudes except in the immediate threshold 
region. Naive extrapolation to the resonance region yields cross sections that are much 
smaller than the available experimental cross sections [[16], O. 



4 Next-to-leading order 

At 0(p A ), the amplitudes consist of the usual two parts: the first one from one- loop 
diagrams with vertices of the lowest-order Lagrangian ( |3.i| ) and a second one from tree- 
level diagrams with exactly one vertex of the chiral Lagrangian of 0(p 4 ): 

^4(4) = ^(4)loop + ^(4)tree • (4-1) 

The loop amplitudes are calculated from diagrams of the general form displayed in 
Fig. where a virtual photon must be appended wherever possible. For the loop amplitude 
we have used a compact representation of the one-loop generating functional for chiral 
SU(2) with at most three propagators [|18| (of O(0 6 ) in the notation of Refs. [0, [|). We 
have checked the result in the limit of a real photon (q 2 = 0) by comparing with the 
general formulas for radiative four-meson amplitudes [fL9~| . Even the reduced amplitude 
AL oop is quite lengthy and we refer to Ref. [18j for the explicit form. Our excuse for 



not reproducing it here is that that the one-loop amplitude will play a relatively minor 
role for the cross sections in the experimentally accessible region that we consider in this 
paper (0.65 < E(GeV) < 1.05). 

The relevant part of the chiral SU(2) Lagrangian of 0(p A ) is given by [Q] (in the 
notation of, e.g., Ref. J|]) 

A = l j(u^} 2 + l j(u fl u l/ }(uV} + ^( X+ } 2 

+ %»X») + ^UTK,u u ]) . (4.2) 
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The low-energy constants (LECs) li,l 2 appear in the (radiative) hit amplitudes, / 3 ,/ 4 
enter through renormalization of the pion mass and decay constant and Iq governs the 
pion charge radius term in the expansion of the pion form factor. The diagrams are the 
same as in Fig. [I] except that exactly one vertex from (|4.2| ) must be inserted, with at 
most one other vertex from the lowest-order Lagrangian ( |3.1| ). The result expressed in 
terms of the reduced amplitude A^ tiee is as follows: 

F*Afa tr J Pl ,p 2 ,p 3 ,p 4 ) = 2l 2 {v-tM (4.3) 

+ hl^s 2 - AsM 2 + AM A ) + |(s 2 - 2ht 2 + 2t\ - 8hu + Av 2 - {q 2 - 2t 3 ) 2 ) 



+2l 3 M A + 2h(sM 2 - M A ) + kq 2 (M 2 - s)} 



P3 



2t 3 -q 
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We have included in A^ tme the chiral logs from the loop diagrams. The quantities k, the 
amplitude (|4.3|) and the loop amplitude A^ loop are then separately scale independent: 

1 M 2 1 M 2 1 M 2 

r 4 = ; . w _J_ ln ^, r^^)^,,^. ( 4.4) 

For the numerical analysis, we use the following values for the LECs that correspond 



to the one- loop analysis in Ref. [20 



Zi = -2.0 x 10~ 3 , l 2 = 1.1 x 10~ 2 , l 3 = -4.6 x 10" 3 , 

k = 2.8 x 10" 2 , f 6 = -1.7 x 10- 2 . (4.5) 

The 0(p A ) cross sections constructed from the amplitude 

A" = Afa + Af 4)loop + Af 4)tree (4.6) 

are shown as dotted curves in Figs. [?],§ for the energy range 0.65 < E(GeV) < 1.05. Com- 
parison with the available data for the 27r + 27r~ channel [jTJjl indicates that the theoretical 
cross sections are still too small. The reason is clear: the amplitudes of 0(p A ) contain only 
the low-energy remainders of resonance exchange. We have to include meson resonance 
exchange explicitly if we want to extrapolate the chiral amplitudes to the 1 GeV region. 



5 Resonance amplitudes generated at 0(p 4 ) 

The renormalized LECs Z[(/x) as well as their SU(3) counterparts are known to be domi- 
nated by meson resonance exchange at typical scales /i ~ M p . The tree-level amplitude 
( |4.3|) of 0{p A ) therefore specifies how to extend the amplitude of 0{p A ) into the resonance 
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region. In the SU(3) notation, the relevant part of the resonance Lagrangian is given by 
C[V(l-),A(l ++ ),S(0 ++ )] = C kin [V,A,S,S 1 ] (5.1) 

+c d (Su^) + c m (Sx+) + CdS^u^) + c m S 1 (x+) ■ 

The octets of vector and axial-vector mesons V(l ), A(l ++ ) are described by antisym- 
metric tensor fields V^, (see App. |B]). S,Si are the scalar octet and singlet fields, 
respectively. Resonance exchange dominance of the LECs at a scale fi = M p amounts to 
the following relations: 

r 2 ( 7$ r 2 \ 



M 2 \M 2 6M 2 fo 



V 2 {M P ) 



Lry 



(~2 ~ ~ 2 



X M 2 6M%J 

H(M P ) - ~ F -^f. (5.2) 

We have omitted the small contributions from kaon and eta loops @] . The axial coupling 
Fa does not enter at this order but will be needed in the following section. At the 577(2) 
level, there is of course no distinction between the SU(2) singlet in S and the SU(3) 
singlet Si. We associate the singlet field in S with the f and the SU(3) singlet with the 
putative a meson. The overall contribution from scalar exchange turns out to be very 
small so that the issue of scalar mixing with or without glueballs [21] has no impact on 
our amplitudes in practice. 

We use M p = 0.775 GeV and the following values for the vector couplings F v , Gy- 

F v = 0.14 GeV , G v = 0.066 GeV . (5.3) 

Gy is obtained from the width T(p — > 7r7r)= 0.15 GeV. The chosen value for Fy is the 
mean value of two possible determinations from T(p° — > e + e~) and from the pion charge 
radius, respectively ||. These values compare well with the theoretically favoured values 

F v = V2F W = 0.13 GeV , G v = F w /V2 = 0.065 GeV . (5.4) 

In the scalar sector, we use || 

c d = 0.032 GeV , c m = 0.042 GeV , q = c,/V3 (i = d,m) , (5.5) 
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Figure 3: Resonance exchange diagrams contributing to the amplitudes of 0(p A ). The 
double lines denote p, a and / mesons. 



with the latter nonet relation holding in the large- N c limit. 



In the tree-level amplitude ( |4.3[ ) of 0(p ) the renormalized LECs /[(M p ) are now set 
to zero and only the chiral logs of Eq. ( fOj) are kept. Instead, the explicit resonance 
exchange diagrams in Fig. ^| are calculated giving rise to amplitudes A^ and A$, thereby 
matching the 0(p 4 ) amplitude on to the resonance region: 

AG 2 fPsPi ■ P2 ~ P1P2 ■ Ps , 2p^[pi • p 4 (p 2 ■ p 3 - t 2 ) - Pi • p 2 (P3 ■ Pa - k)] } 



D p [(pi+p 3 y 

+F v G v (p%t l -p>£t 3 )^ 

F v G v g\s-Ml) „ 
+ ZW)(2t 3 -g 2 ) Ps 



(2t 3 -qi)D p [(p 2 +p 4 )*} 
1 1 



D p (qi) D p [( Pl +p 3 ) 



(5.6) 



S=fo,c 

with scalar couplings 



2[{s-2M$)c* + 2Ml c 
3(2i 3 - q 2 )D s (s) 



T2„S 12 



" J 



Jo 



y/E't 



d, m) 



and propagators with energy dependent widths 



D P (t) 

r,(t) 
r*(t) 



Mp-t- iMpT P (t) 



M„t 



^ f(l-4M*/f) 1/a 



4M 2 /f) 3/2 #(t - 4M 2 ; 



Mf(l 



4M2/M|)V 
9 



(5.7) 

(5.8) 

(5.9) 
(5.10) 

(5.11) 



In the numerical analysis, we take 



M fo = 0.98 GeV , T /o = 0.05 GeV , 

M a = 0.6 GeV , r CT = 0.6 GeV . (5.12) 

Finally, we recall that FyGy and CdC m are both positive f|. 

At this point, the amplitude A^(pi,p2,P3,P4) takes the following form: 

N 1 = + A (4)loop + A^ tiee + A^ + A^ s (5.13) 

where the amplitude A^ tiee contains only the chiral logs in ( [4.4| ) for p = M p . The 
renormalized LECs /[(M p ) have been traded for the explicit resonance exchange ampli- 
tudes A^,Ag. The resulting cross sections are more realistic than the strictly 0(p 4 ) cross 
sections from the amplitude (£OJ), but they are still too small in comparison with the 



available data around 1 GeV [[Lq, 17 



6 Beyond 0(p z 

There must be additional ingredients in the amplitudes that make important contribu- 
tions to the cross sections already for energies below 1 GeV. To the extent that the LECs 
/[(M p ) are known to be dominated by p (and scalar) exchange as given in Eqs. ( |5.2| ). 
those additional amplitudes must vanishQ to 0(p 4 ). We will try to locate the dominant 
contributions that appear first at 0(p 6 ) but in contrast to the previous section we cannot 
claim completeness here. At this order, also diagrams with more than one meson reso- 
nance contribute. Even for single resonance exchange, a complete analysis of 0(p 6 ) is not 
available at present. 

However, we may turn to existing phenomenological treatments |9|, [H| [TlJ for guid- 



ance. In addition to the obvious p (and the less important scalar) exchange, the data 
16, [lTj clearly indicate the presence of u and a\ exchange. 



6.1 Double p exchange 

The resonance Lagrangian (|5.1| ) also generates amplitudes starting at 0(p 6 ) with two 
p mesons exchanged. The corresponding diagrams are displayed in Fig. |j. The diagram 
where the virtual photon couples to p + p~ is actually required by gauge invariance because 
the (charged) vector mesons are dynamical fields here. The diagrams of Fig. |3| produce a 
gauge invariant amplitude in the strict 0(p 4 ) limit only where the resonance propagators 
shrink to points. Although of different chiral order, the amplitudes of Figs. |3| and |] must 
be added for a meaningful amplitude in the resonance region. 

2 Small additional contributions to the LECs of 0(p 4 ) are possible and even expected, e.g. from p 1 
exchange. 
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Figure 4: Double p exchange diagrams generated by the Lagrangian (|5.1|) . 

All couplings needed for the diagrams of Fig. |4] have already been defined. The (re- 
duced) double p exchange amplitude has the explicit form 



F n A pp(Pl,P2,P3,P4) = 

v21 {p 2 (M 2 + Pi ■ P:i)(Pi - Ps) ■ Pa 



4G 2 V 



D P [(pi+P3) 2 ]D p [(p 2 +p 4 ) 

+p%{Ml + p 2 ■ p 4 )pi • (P2 - Pa)} + (6.1) 

DM 2 )D p [{p\ + Pa?} ^ q2(Pl " P3) ' P4 + ^ ' ^ ~ P0 

' [(Pi - Psf (*4 - t 2 ) - (P2 - P4)^(t 3 - tl)](M^ + p 2 • p 4 )} + 

8F V GW 



F2D p (qZ)D p [(p 2 +p 4 )*](2t 3 - q 2) 



Ps {Pi ■ Pa(P2 ■ Ps - h) - Pi ■ P 2 (P3 ■ Pa - U)} 



Putting together the lowest-order amplitude A'jU in (|3.3| ) and the single and double 
p exchange amplitudes Ap 1 , A^ p in (|5.6j), (|6.1|), one observes that some terms can be 
combined as the leading terms in the low-energy expansion of the pion form factor in 
view of the relation F V G V — F% 0: 

F (a 2 ) - 1 + FvGvq2 + J_r (a 2 ) + ~ Mp2 (6 2) 

Mg j " + f 2 /w) + m; m ] + ■ ■ ■ - dm 2 ) ■ { } 

We have included the leading-order absorptive part iT p (q 2 )/M p that is contained in the 
one-loop amplitude A^ loop in one case and is of higher order in the other case. Replacing 
the expansion terms on the left-hand side of ( |6.2j ) by the usual p dominance form (the 
right-hand side of Eq. (|6.2j)) is equivalent to the order considered. Of course, the partial 
resummation yields a phenomenologically much more realistic amplitude. A similar re- 
summation applies to the scalar exchange amplitude. We therefore express the combined 
p and scalar exchange amplitude in the following form 

(Afa + A$) M 2 /D p (q 2 ) + + A» p . (6.3) 
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The modified single and double p exchange amplitudes are now given by 



AG 



2 \vlVl ■ P2 - PlP2 ■ P3 , 2 Ps [Pi • P*(P2 ■ P3 - h) ~ Pi ■ P2(P3 ' P i ~ U)}M* ' 

) 



\ D p [(p 1+P3 y] (2t 3 -qi)D P i(p 2 +p 4 y]D, 



+F v G v (p%t 1 - I f 1 t 3 ) (-1^ - - 1 I (6.4) 



FnA p l p (p 1 ,p 2 ,P3iP4) = 

yi] [p^{Ml + pi ■ p 3 )(pi - p 3 ) • p 4 



4G 2 , 



+Ps( M ^ + P2 ■ P4)pi ■ (P2 ~ Pa)} + 

+ [(Pi - p 3 )"(*4 - h) - (p 2 - p 4 f (t 3 - ti)](M2 + p 2 • p 4 )} . (6.5) 

The leading-order absorptive part of the pion form factor must be subtracted from the 
one-loop amplitude yielding a modified loop amplitude A^ loop . 

6.2 u; exchange: vector meson dominance 

Vector meson dominance (VMD) for u decays postulates the dominant role of an upn 
coupling p2| ^3fl . We write the corresponding Lagrangian (unique to lowest order in 
derivatives) as 

C(ujpn) = g U pn£^ P aOJ IJ 'd u 7f ■ pP a . (6.6) 

In this case, it is more convenient to describe the u in terms of a vector field (see 
App. i). 

The decay u — > p°7r° — > 7r°7 proceeds with a rate 

,2 T?2 



r(cu - 7r° 7 ) = (m 2 , - M 2 ) 3 . (6.7) 



6M 3 M 



The measured partial width |24| corresponds to \g Ufm \ = 5.0. On the other hand, the 
dominant decay chain uj — > pit — > 3tt leads to \g wpn \ = 5.7. A small direct uo — > 3tt 
amplitude is allowed but VMD clearly accounts for the dominant features of both decays. 

For e + e _ — > 47r, the relevant uo exchange diagram is shown in Fig. It gives rise to 
a (reduced) amplitude 



A£(pi,P2,P3,P4) = „ 9 Fi 8F ^ V n L °' m ^t{-P2*4P1-P3+P3(*4P1"P2-*2P1-P4)} 



2 

F^D p ^)D u [{q-p 1 f] 
{D; l [(p 2+P3 ) 2 ] + D- p l [(p 2 +p,f] + D- p l \(p 3 +p 4 ) 2 ]} . (6.8) 
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v/VWW 




Figure 5: VMD diagram for ui exchange. 
In view of the small value = 8.44 MeV [24] we employ an energy independent width 



in the propagator function Du[(q — pi) 2 }. The amplitude ( |6.8|) completely dominates 
the cross section for e + e~ — > 27r°7r + 7r~ around 1 GeV in accordance with experimental 
findings |17|, [TT[] . In order to appreciate the size of this amplitude of 0(p 6 ), we compare 
it to a typical p exchange amplitude of 0{p A ) as given in (|5.6|). By naive chiral counting, 
the dimensionless quantity c w defined by 



24a 2 F 2 



(4^)2 MIMl 

would be expected to be of 0(1). With \g Ufm \ = 5.7 one finds instead \c w \ = 24, quite a 
drastic deviation from naive chiral counting. The sign of the u exchange amplitude (|6.8j ) 
is fixed by the positive sign of FyGy [§, [|. Of course, the corresponding amplitude due 
to cf) exchange is completely negligible. 

6.3 ai exchange 

Although u exchange dominates the cross section for the 27r°7r + 7r~ final state it does not 
contribute to the other channel 2ti + 2ti^ . Here a\ exchange will play an important role. 
We follow the usual VMD assumption that the dominant decay mode a\ — > 3n proceeds 
via an intermediate p. 

Contrary to upir, there are several possible chiral couplings for the aipn vertex. The 
ideal place to analyse this vertex is the process r — > z/ r 37r and such an analysis is under 
way ||25|| . In the tensor field formalism, there are altogether five a priori independent aipn 



couplings of lowest possible chiral order [^j . Two of them give the same amplitudes in 
our case and another one is proportional to M 2 and therefore vanishes in the chiral limit 



25]. We will restrict ourselves here to the remaining terms that boil down to the following 
Lagrangian for the charged a\ fields (the neutral a\ cannot be exchanged in the diagrams 
of Fig. 0): 



M ai 

+ ^a + r (d,pld\~ - <V^«9V) + h.c. , (6.10) 
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with dimensionless couplings c 2 ,c 3 ,c 4 . 

The analysis of r — > z/ r 37r that should determine or at least relate the constants c, is 
not yet available |[25|| . In order to proceed, we make the simplifying assumption that the 
couplings are all equal: 

c 2 = c 3 = c 4 . (6.11) 

From the decay width T(ai — > p7r — > 37r), accounting for the finite p width, we find 

|c 2 | = 319 [r(ai -» 3tt)/0.5 GeV] 1/2 . (6.12) 

The surprisingly large value of | ca | is due to the fact that the (off-shell) a\ — > pit — > 37T 
vertex function vanishes in the chiral limit for the choice (|6.11|) , as long as the pions 
are on-shell. Although this property is certainly not required by chiral symmetry it has 
interesting implications for the high-energy behaviour of the r — > z/ T 37r amplitude . 
In this case, the neglected coupling proportional to M% should be added in a complete 
analysis of the a\pit vertex. We have also investigated other choices for the couplings q: 
the resulting cross sections are always smaller than for the choice Q6.ll 




x/WW 




Figure 6: a\ exchange diagrams. 

As shown in Fig. || there are two a\ exchange diagrams that must be taken into 
account here. The first one has a direct aiifj coupling Fa defined in the resonance La- 
grangian (|5.1|) . We take the theoretically favoured value |j| — F n for this coupling. 
In addition, also the a^pit couplings in ( |6.1U| ) contribute to the radiative decay a\ — > ify. 
For the choice ( |6.11|) we obtain an effective coupling 



Fa ~ Fa ~ M ai M* ' (6 - 13) 

with the two terms approximately equal in magnitude. Requiring constructive interference 
(sgn(c 2 i 7 AFyF 7r ) < 0), the resulting partial width T(ai — > ttj) is larger than the PDG 
value of 640 keV (based on a single experiment) by about a factor 2.5 for r(ai — > 3n) = 
0.5 GeV. 

The (reduced) amplitude ^^(^1,^2,^3,^4) from the two diagrams in Fig. || is in an 
obvious notation given by 

M l D ai [{q ~ Pa) 2 \ D p (Pi + 
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V R 



M ai D p (q*) 

-T-7T 



9fip9v<r(M ai k ) + g^pk y k a g^kykp (/i 
F 2 M ( 

+ C4P2 • (pi +Ps)P3Pl} 



— {c 2 (pi • - P2 ■ P3P2P1) + ca(pi • Pa + Ml)p p 2 (pl - pf) 

ai 



The relative signs in this amplitude are determined by the choice (|6.11|) , the constructive 
interference in ( |3.13| ) and by FyGy > 0, ||. 

The energy dependence of the a\ width T ai (t) has also a considerable numerical im- 
pact. Awaiting the results of Ref. [p5l, we assume for the present analysis the functional 
form (all numbers are to be understood in appropriate units of GeV) suggested by Kiihn 
and Santamaria 



T ai (t) = r ai g(t)/g(M 2 ai ] 
g(t) = (1.623 1 + 10.38 - 9.23/t + 0.65/t 2 ) 6[t - (M p + M^) 2 } (6.15) 
+ 4.1(t - 9M 2 ) 3 [1 - 3.3(t - 9Af 2 ) + 5.8(i - 9M 2 ) 2 ] &(t - 9M 2 )Q[(M p + M n ) 2 - t] . 




Figure 7: Comparison of data JTB| (left figure) and predictions (right figure, see text) for 
the cross section a(e + e" -> 2n + 2n~) for 0.65 < £(GeV) < 1.05. 

Putting all contributions together, we arrive at our final amplitude 

^nal = K 2 ) + A%} M 2 p l DM 2 ) + ^4)loop + ^tree + K + K P +K + K ■ (6-16) 
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0.6 0.65 0.7 0.75 O.S 0.85 0.9 0.95 1 1 .05 



energy (GeV) 

Figure 8: Theoretical predictions for the cross section er(e + e~ — > 27r°7r + 7r~) for 0.65 
< E(GeV) < 1.05. 

This amplitude has the correct low-energy behaviour to 0(p 4 ) and is expected to contain 
the relevant ingredients for an extrapolation to the 1 GeV region. 



7 Cross sections and decay rates 

In Fig. [7|we compare our results for the cross section a(e + e~ — > 2tt + 2tt~) with available 



data taken from Ref. [16|. The dotted curve corresponds to the strictly 0(p A ) amplitude 
( Of ) and the full curve is the cross section for the final amplitude ( |6.16| ). Whereas the 
dotted curve is definitely too low the full curve agrees well with experimental data up 
to 1 GeV. The more pronounced rise of the full curve is mainly due to a\ exchange. The 
p exchange amplitude generated at 0{j> 4 ) is also important, to a lesser extent also the 
lowest-order amplitude with resummed pion form factor. Loops and chiral logs are much 
less important. Finally, scalar exchange contributes very little to the cross section and 
double p exchange does not contribute at all in this channel. 

At low energies where our amplitude should be most reliable the predicted cross 
section is below the shaded band in Fig. [7|. That band [16| describes an extrapolation 



from data at higher energies [17] with a non-chiral resonance model. 

The analogous theoretical results for a(e + e~ — > 27r°7r + 7r~) are shown in Fig. |8|. Al- 
though there are so far no data in the region below 1 GeV the theoretical cross section 
(full curve) connects well with the data starting at 1 GeV |T7|]. For this channel, the much 
steeper rise compared to the 2-n + 2n~ mode is almost exclusively due to to exchange. a\ 
exchange, even though less important here, interferes constructively with u exchange in 
the 1 GeV region. All other contributions are very small there. 

Near the p pole, the amplitudes are of course completely dominated by the resonant 
parts containing the propagator function 7J^ 1 (g 2 ). The cross sections at E = M p therefore 
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determine the branching ratios for p — > An according to the general formula 

For the channel 2tt + 2tt~, the relevant contributions are the lowest-order amplitude with 
pion form factor, single p and a\ exchange. There is a destructive interference between the 
modified lowest-order amplitude on the one hand and the two single-resonance exchange 
amplitudes at E = M p . Except for a\ exchange, this interference is dictated by the QCD 
structure of 0{p A ). Although the a\ amplitude depends on our assumption ( |6.11| ) for the 
a\pn vertex the relative sign to the other two amplitudes is also fixed. The situation 
is a little different in the 2n°n + n~ channel because of the additional uo exchange. In 
this channel, the interference pattern is: lowest-order amplitude with pion form factor 
+ di exchange - single p exchange - u exchange. Although the couplings involved are 
relatively well known we assign a 40 % error to the calculated branching ratios in view 
of the destructive interferences: 

BR(p° -> 2n + 2n~) = (6.7 ± 2.7) x 
BR(p° -> 2tt°ti + tt-) = (5.0 ± 2.0) x 10~ 6 . (7.2) 



The result for the 2n + 2n mode agrees within errors with the experimental value pi | 



(extracted from the cross section in Fig. [7] at E = M p ) although our mean value is almost 
a factor three smaller. For the 2tt q tt + tt~ channel there is only an experimental upper limit 
f24}j that is compatible with ( [7.21) . There is a wide range of model predictions for the An 



decay modes of the p as reviewed in Ref. |27 



We have plotted the cross sections only for energies < 1 GeV because our amplitudes 
do not have the correct high-energy behaviour. This manifests itself already in the 1 to 
2 GeV region where our cross sections exceed the experimental cross sections |T7| 



Scaling all four-momenta in the same way and requiring that a(e + e~ — > An) decreases 
at least as fast as (most likely faster than) 1/E 2 at large energies to satisfy the asymptotic 
QCD constraint, one finds that the basic current matrix element J ,M (pi,P2,P3,P4) in (|2~6l ) 
must vanish at large energies at least as 1/E. This criterion is not even met by the lowest- 
order matrix element ( [3.4|) although it is satisfied by the modified lowest-order amplitude 
in (|6.3| ) due to the pion form factor. 



In addition to resummations of parts of the amplitude, additional higher-mass states 
must be included in order to access the region up to 2 GeV and to satisfy the high- 
energy constraints of QCD. The Particle Data Group lists many such resonances with 
the appropriate quantum numbers, e.g., p(1450), p(1700), /o(1370), /o(1500), /o(1710), 
a;(1420), a;(i650) and states with higher spins. In the spirit of duality, all those states 
are expected to conspire to produce the right asymptotic behaviour of the amplitudes at 
high energies. 
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8 Conclusions 



We have performed the first calculation of the processes e + e~ — *> Air and r — > z/ T 47r with 
the correct structure to 0(p 4 ) in the low-energy expansion of the standard model. In 
addition to the proper low-energy structure, CHPT automatically produces amplitudes 
with all relevant symmetries of the standard model, in our case (spontaneously and softly 
broken) chiral symmetry, gauge invariance, Bose symmetry and C invariance. 

Although the chiral amplitude to 0(p A ) is only valid close to threshold it contains 
information how to extrapolate to the resonance region. This information on p and scalar 
exchange is however not sufficient to describe the e + e~ cross sections up to energies of 
around 1 GeV. We have therefore included as additional contributions to, ai and double 
p exchange that first show up at 0(p 6 ). All necessary couplings were determined from 
the decay widths of the various resonances involved. 

The predicted cross sections for E < 1 GeV and the branching ratios BR(p° — > Att) 
are in good agreement with available data. Our amplitudes do not have an acceptable 
high-energy behaviour so that additional ingredients are needed (such as higher-mass 
resonance exchange) to make predictions in the phenomenologically interesting region up 
to 2 GeV. 

In the isospin limit, the r decay amplitudes can be calculated from the annihilation 
amplitude for the channel 27r°7r + 7r~ fl3|| . The comparison with r decay data will be 
postponed until the proper high-energy behaviour has been implemented. 
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A Isospin relations 



From the isospin relations [13| (2.6),. . . ,( [2.9|) for the four possible final states it is evident 



that the amplitude for the 27r°7r + 7r~ channel is sufficient to obtain the remaining three 
amplitudes. From the observation that u exchange cannot contribute to the 27r + 27r~, 
37r°7r _ modes one finds immediately that none of those latter modes is sufficient to cal- 
culate the remaining ones. 

The only nontrivial question?] concerns the channel 27r _ 7r + 7r°. At first sight, one would 
expect that from the sum of two current matrix elements in ( |2.8| ) one cannot determine 
J tJ- (pi,p 2 ,P3,P4) itself. However, one has to take into account the symmetry relations for 
J^ipiiPiiPsiPij- It is the purpose of this appendix to show explicitly that knowledge 
of the amplitude for the 27r~7r + 7r° mode is also sufficient to determine all four matrix 
elements in the isospin limit. 

For this purpose, we write the general decomposition of J^(pi,P2,P3,Pi) as 

J^(Pl,P2,P3,P4) = PlB(pi,P2,P3,P4) +P2 B (P2,Pl,P3,P4) 

+ P3 C (PliP2,P3iP4) ~ P4C(p U p 2 ,P4,P3) (A.l) 

in terms of two invariant amplitudes that satisfy the constraints 

B(pi,P2,P3,PA) = -B(pi,p 2 ,P4,P3) 

C(pi,P2,P3,Pi) = C(p 2 ,Pl,P3,P4) (A.2) 

due to charge conjugation invariance and Bose symmetry. Gauge invariance leads to a 
further relation between B and C but we do not need this relation here. 
The isospin relation ( |2.8|) can now be written as 



(7T-7r-7r + 7r°|V^(0)|0>/V2 = P^D(p + ,p u p 2 ,p ) +p^F(p + , Pll p 2 ,p ) 

+ P2F(p+,p 2 ,Pi,Po) ~PoG(P+,Pi,P2,Po) (A.3) 

with 

D(p+,Pi,P2,Pa) = B(p+,pi,p 2 ,po)+B(p + ,p 2 ,pi,po) 
F(p+,pi,P2,Po) = B(pi,p + ,p2,po) + C(p + ,p 2 ,Pi,Po) (A.4) 
G(p+,Pi,P 2 ,Po) = C(p+,PuPo,P2) + C(p + ,p 2 ,p ,p 1 ) . 

With the symmetry relations ( |A.2| ) one easily verifies 

2B(p 1 ,p 2 ,p 3 ,p 4 ) = D(px,P3,P2,Pi)+F(p2,Pl,P3,P4)-F(p3,Pl,P2,Pt) 

2C(px,p 2 ,p 3 ,p A ) = -D(p 3 ,p 2 ,pi,p A ) + F(p 1 ,p 3 ,p 2 ,p A ) + F(p 2} p 3l pi, p 4 ) . (A.5) 

Therefore, the amplitude for the 27r°7r + 7r _ channel can be obtained from the 27r~7r + 7r° 
amplitude. Consequently, knowledge of the 27r~7r + 7r° mode is enough to determine the 
other three amplitudes in the isospin limit. 

3 We thank Hans Kiihn for raising this question. 
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B Vector and axial- vector mesons 



Spin-l mesons can be described either by the more conventional vector (or axial- vector) 
fields or by antisymmetric tensor fields V^ u . For matching resonance exchange ampli- 
tudes with standard CHPT amplitudes, the choice of fields is a matter of convenience. 
The tensor field formalism has the advantage of producing immediately the correct LECs 
of 0(p 4 ) 0, ||. On the other hand, single resonance exchange that contributes first at 
0(p 6 ) such as us exchange is better described by vector fields |5|, P% . The two descriptions 
are equivalent but the transformation from one formalism to the other involves the in- 
troduction of explicit local amplitudes. Those local terms may be avoided by employing 
the proper formalism. 

We recall first the usual normalization of a conventional massive vector field for a 
vector meson of mass M, with polarization vector e^p), and the associated propagator: 

(0K(0)\V,p) = e» (B.l) 

(2n)*(AP-k*-ie) ^ " KKI ^ ' ^ 

The same spin-l particle can also be described by an antisymmetric tensor field V^. 
The corresponding one-particle matrix element and the propagator are given by [|J 

(0\V^(0)\V,p) = iM~ 1 {p fl e u (p) —PuBfi{p)} (B.3) 

(2ir)*(M* - V - ie) W" 9 ^ ~ ^ 

+ gwKK - QfiaKkp - (fx «-> v)\ . (B.4) 

In many cases, the tensor field propagator can be simplified. Whenever the p meson 
couples either directly to the (virtual) photon or to two pions the transverse part of (p.4| ) 
does not contribute [|7| and the p propagator may be replaced by 

(0\T{ P (0)}|0> = i 

(27r) 4 (M 2 — k 2 — is) ^^ va 9ftcr9up] ■ 

(B.5) 

This happens to be the case for all diagrams considered involving p exchange. The sim- 
plification does not apply for the a± propagator in the diagrams of Fig. |6|. 
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C Numerical input 



In this appendix we collect the numerical values of masses and coupling constants that 
we have used for the calculation of cross sections. 

Chiral LECs 

F n = 0.0924 GeV \ = -2.0 xlO" 3 [ 2 = 1.1 xl0~ 2 

/ 3 =-4.6xl0- 3 / 4 = 2.8xl0- 2 l 6 = -1.7xl0- 2 

Vector mesons 

M p = 0.775 GeV F v = 0.14 GeV G v = 0.066 GeV 

M w = 0.783 GeV Y u = 0.00844 GeV g Ufm = 5.7 

Axial-vector meson 

M ai = 1.23 GeV T ai = 0.5 GeV 

Fa = F n c 2 = c 3 = c 4 = 319 

Scalar mesons 

M fo = 0.98 GeV T /o = 0.05 GeV 

M a = 0.6 GeV T CT = 0.6 GeV 

c d — 0.032 GeV c m = 0.042 GeV q = a/y/3 (i = d,m) 
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